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Q,, Maxwell's equations in curved space-time are invariant under electromagnetic duality transfer- 
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mations. We exploit this property to constraint the design parameters of metamaterials used fer 
transfermations optics. We show that a general transfermation must be implemented using a dual- 



jy^' . symmetric metamaterial. This can be accomplished constraining the polarisability tensors of their 

O 

individual constituents, i.e. the meta atoms. We obtain explicit expressions for these constraints. 

O ' We also show that the spatial part of the coordinate transformation depends only on the electric- 

O ■ electric tensorial coupling of the polarisability tensor, while the spatio-temporal part depends only 

J>-j, on the electric-magnetic tensorial coupling. In our derivations, we find that two dipoles located at 

Qh| the same point, one electric (d) and one magnetic (m), are needed to produce a total field with 

well defined helicity equal to +1 or -1, and that they must be related as d = ^m or d = — ^m, 
> ■ 
Q"^ . respectively. 
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Transformation optics offers us a path to the design of invisibihty cloaks, perfect lenses 
and any other device whose action on the electromagnetic field can be casted as a coordinate 
transformation in space-time l|-l3|. Transformation optics is based on the exact correspon- 
dence between Maxwell's equations in an empty region of curved space-time and inside a 
material medium in a flat space-time background. All the information about the curved 
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space-time geometry is contained in the constitutive relations of the material j^. The de- 
sired transformation specifies a space-time metric which at its turn specifies the constitutive 
relations of the material. Refer to 5| for a detailed treatise in transformation optics. 

Such formidable step in our ability to manipulate light comes with a correspondingly 
steep increase in the tunability requirements of material constitutive relations. Nature does 
not provide us with nearly enough flexibility in this aspect. We must synthesize artificial 
materials: Electromagnetic metamaterials. Transformation media are typically implemented 
by means of electromagnetically small inclusions inside an homogeneous and isotropic di- 
electric. These inclusions are sometimes referred to as meta atoms. The idea is to obtain the 
required constitutive relations from the collective response of the meta atoms. Currently, 
though, there is no systematic design methodology to go from the constitutive relations to 
the actual implementation of the metamaterial. In general, this is a highly complex task, 
partly because of the huge number of degrees of freedom of a general metamaterial, which 
include the electromagnetic response of the meta atoms and their three dimensional spatial 
arrangement. Reducing the number of degrees of freedom while maintaining the ability to 
implement general coordinate transformations is desirable. 

In this paper we use a non-geometrical symmetry of Maxwell's equations, electromagnetic 
duality, to constrain the individual response of the meta atoms. There is a deep connection 
between transformation optics and duality symmetry. Almost two decades ago, I. Bialynicki- 
Birula realised that that the two helicity components of light do not mix in a gravitational 
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Since helicity, as an operator, is the generator of duality transformations in the 
same same sense that angular momentum is the generator of rotations, helicity preservation 
is equivalent to invariance under duality transformations. It follows that Maxwell's equa- 
tions in a general space-time geometry are invariant under duality transformations. We can 
therefore obtain a constraint without sacrificing generality: a metamaterial that implements 
a transformation medium should preserve the helicity of light, or in other words, possess 
duality invariance. In this paper, we will assume that the meta atoms in the metamaterial 



are small compared to the electromagnetic field wavelenghts so that they can be treated in 
the dipolar approximation. We will study the restrictions that the requirement of duality 
invariance imposes on the polarisability tensors of the meta atoms. Electromagnetic dual- 
ity is a transformation that mixes electric and magnetic fields by means of a real angle 9. 
Assuming space and time dependent fields (r,t): 



E^Ee = Ecos^-ZoHsin^, 

(1) 
ZqH -^ ZoUe = E sin + ZqU cos 6, 



where Zq = \/Jio/eo and (eo,/io) are the vacuum permittivity and permeability con- 
stants. In vacuum, ([T]) is a symmetry of Maxwell's equations: If the electromagnetic 
field (E(r, t), H(r, t)) is a solution of the free space Maxwell's equations, then the field 
r, t), H6)(r, t)) is also a solution for all 6. Some decades ago. Calkin |8|, and then others 



10| ) showed that helicity was the conserved quantity related to such symmetry. The 



helicity operator is defined H, chap. 8.4.1] as the projection of the total angular momen- 
tum J onto the linear momentum direction, i.e. A = J ■ P/|P|. In the case of photons [l^, 
chap. 2.5], helicity can take the values ±1, which completely describe the internal degrees 
of freedom (polarisation) of the photon and the electromagnetic field. In the same way that 
linear momentum generates translations and angular momentum generates rotations, the 
helicity operator generates the duality transformation in ([T]). 

In general, the presence of matter breaks the symmetry of the equations: a solution 
(E(r, t), H(r, t)) does not generate a new solution when transformed as in ([T]). As a conse- 
quence, the interaction of light with matter generally changes the helicity of light. Neverthe- 
less, the electromagnetic duality symmetry can be restored for the source-free macroscopic 
Maxwell's equations in material systems characterized by scalar permittivities and perme- 



abilities 



Ci. 



We shall see that duality can also be restored in a macroscopic bianisotropic 
and inhomogeneous medium and in the dipolar approximation. Both of these cases are 
relevant for transformation devices made with metamaterials. 

In this paper, we will use the Riemann-Silberstein representation of electromagnetic fields 
|6|. This formulism is very well suited to treat problems involving duality and helicity: As 
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discussed in [7|, §2.1], assuming that the energy of the field is positive, the combinations 

F±(r,t) = ^(ZoD(r,t)±^B(r,t)) 

f (2) 

G±(r,t) = -^(E(r,t)±zZoH(r,t)). 

separate the positive and negative hehcity components. That is, F_|_(r,t) and G+(r, t) have 
a well defined helicity equal to 1 and F_(r,t) and G_(r, t) have a well defined helicity 
equal to minus 1. Indeed, it is easily seen that under the duality transformation in ([1]), 



G^(r,t) -^ exp 
transformation 
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^0)G±{r,t) and that F^(r,t) — )■ exp(=i=26')F-|-(r, t) under the companion 
chap. 6.11] (the (r,t) dependence now implicit): 

ZqD -^ ZoDe = ZqD cos ^ - B sin 6, 

(3) 
B-^Bg = ZoDsin^ + Bcos^. 



It follows that F-|-(r, t) and G-|-(r, t) must be eigenvectors of the generator of the contin- 
uous transformations ([3]) and ([1]) with eigenvalues equal to ±1. Under the assumption of 
positive energies, such generator is the helicity operator, as we now show. In the free space 
Maxwell's curl equations written as Tj: 

^afF±(r,t) = ±coVxF±(r,t), (4) 

where Cn = l/^eo/xo- We identify idt with the Hamiltonian "H and Vx as the operator 
XII. 93], that is, the projection of the spin angular momentum S onto the linear 
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momentum P. Then, (jlj) can be written like 

co^F± = ±F±. (5) 

If we now assume that the energy is positive, and since we are talking about massless 
electromagnetic fields, then "H = Co|P|. By additionally noting that J-P = (rxP + S)-P = 
S • P, the helicity operator A = J • P/|P| appears explicitely in ([5]): 

AF±(r,t) = ±F±(r,t). (6) 

In connection with the previous discussion about duality, @ shows that helicity has, for all 
F-|-(r, t), the behavior corresponding to the generator of duality transformations. Finally, it 
is important to note that the assumption of positive energy does not reduce the generality 



of the treatment because, for electromagnetic fields, all the information is contained in the 
positive energy part and duplicated in the negative energy part, or vice versa. 

The same arguments hold for the case of a homogeneous and isotropic medium with 
constitutive relations D(r, t) = eE, B(r,t) = yuH(r,t) by simply replacing Cq and fiQ by e 
and /i, hence Zq by Z and Cq by c. 

To exploit the simple transformation properties of F±(r,t) and G±(r, t) under duality, 
we will use the following two six-component vectors: 



V2 



ZT> + iB 
ZT> - iB 



^-T. 



E + iZU 
E-iZU 



(7) 



In (ITj), and from now on, the (r, t) dependence of the fields will be implicitly assumed but not 
explicitly written most of the time. Both vectors in ([7]) keep the positive helicity component 
on the upper three elements and the negative helicity component on the lower three elements. 
We are now ready to start the derivations. 

For starters, we examine the duality transformation properties of Maxwell's equations in 
a medium with general linear constitutive relations J-" = NQ and derive the conditions on 
A^ for the medium to be dual symmetric. We implicitely assume that A^ can in general be a 
function of space. Then, we focus on the particular case of the linear constitutive relations 
induced by a general space-time geometry and verify that they meet the duality conditions. 

Using the curl equations we can write: 

Vx 

-Vx 

In the Q{J^) basis, the duality transformation of equations ([1]) and ([3]) is simply the 6x6 

matrix: 

exp{-ie)I 

exp(ie)/ 

where / is the 3x3 identity matrix. We now use (jH]) to transform 

Vx 



idtJ' = idtNg 



(8) 



Da 



(9) 



DeidtND^'DeG = Dg 



D^'DeG. 



(10) 



-Vx 

Since Dg commutes with idt and with the block diagonal operator containing the curls, we 

obtain 

Vx 



idtDeND.'g 







-Vx 



(11) 



where Q = DqQ. The necessary and sufficient invariance condition that keeps the form of 
equation (fTTj) in the Q variable the same as the form of equation ([8]) in the Q variable is 
DqNDq^ = N, that is, that A^ and Dg must commute. This happens if and only if A^ is 
block diagonal in 3x3 blocks. If we work back what this block diagonal condition means in 
the constitutive relations among the fields: 



(12) 
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we obtain that a block diagonal A^ forces 

e = 

with which the time evolution of Q reads: 

~iX 



Ai, X 



-7, 



idt 



e + ix 





Vx 





Q = 











-Vx 



(13) 



(14) 



The restrictions in ( TT3l) are the necessary and sufficient conditions for duality invariance 
(helicity preservation) for a general linear inhomogeneous and bianisotropic media. Inci- 
dentally, we recall that the constitutive matrix A^ can be a function of time and space, 
independently of which, if ( TT3|) is met, helicity is preserved. In particular, this applies across 
ihe boundaries of different bianisotropic media, which generalizes the results obtained in 
3, §2.2] for isotropic and inhomogeneous media without magneto-electric couplings and in 



13l | for helicity preservation across boundaries of different isotropic and homogeneous media 



without magneto-electric couplings. 

We now turn our attention to the constitutive relations induced by an empty but curved 
space-time with a metric g^^. They were derived in ^. In the Riemann-Sielberstein repre- 
sentation, they can be shown to be [7]: 



n = -^ (v^^""^ - ^gokS'"'"') Gl (15) 

where n,m = [1, 2, 3], y/—g is the square root of the determinant of —g^y, the inverse metric 
g^^ is such that g^^guk = ^k ~ diag(l, 1, 1, 1), and 5"'^™ is the totally antisymmetric three 
dimensional Levi-Civita symbol. 

For the other helicity components (G_, F_), straightforward algebra leads to: 



( / Tiri'n-m I • nkm\ r~i- 
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(16) 



We can combine (fTS!) and (fT6!) in a matrix form: 



T 



A+ 
A_ 



(17) 



where A!^ = {—^^—gg"'"^ T Wok£"'^"^)/goo- From our previous discussion, the block diagonal 
form of (1171) implies duality invariance of a curved space-time. Since we have not imposed 
any restriction on g^^, helicity preservation is inherent to the structure of any space-time 
metric, and can hence be seen as a necessary condition for any transformation medium, 
furthermore, using results from |3|, we can now arrive at some interesting conclusions. In 
3|, the authors showed that e is related to space-only transformations and x to transfor- 
mations which mix space and time components. In equation ( fT4l) we note that, since the 
curl operator has an opposite sign effect on the two helicity components, e has the same 
effect in the time evolution of both helicity components while x causes a helicity dependent 
transformation. Consequently, we arrive at the conclusion that space-only transformations 
are helicity independent while space time mixing has a helicity dependent effect. 

We will now show how the requirement of duality invariance in transformation optics 
affects the design of metamaterials. Clearly, if all the constituent meta atoms are dual 
symmetric their collective response will also be dual symmetric. In other words, if the 
field scattered by each meta atom preserves the helicity of its exciting field, helicity will be 
preserved by the entire meta medium. There is strong evidence that the converse is also true 
in the general case [16!]: the collective response of a mixture of scatterers preserves helicity 
only if each individual scatterer preserves helicity. To impose duality invariance to the meta 
atoms, we assume that the size of a meta atom allows to model it by a polarisability tensor 
. Under these conditions, we ask what restrictions on the polarisability tensor result from 
the requirement of duality symmetry. From now on, we work on the frequency {u) domain. 
This treatment can better accomodate the frequency dependent response of many practical 
meta atoms. It is understood that the obtained restrictions must be met across the whole 
bandwidth of operation. 

We start with a dipolar scatterer at position r' with polarisability tensor M. The electric 
(d) and magnetic (m) dipoles induced by an incident electromagnetic field can be written 
as the product of the 6x6 real tensor M times the six component column vector formed by 



the incident electric E(r', w) and magnetic H(r', w) fields: 



dicu) 




B{r',Lo) 






= M 




= 


m{uj) 




H(r',u;) 





ttmE ttmH 



E(r',u;) 



(18) 



where M is decomposed into its four 3x3 blocks, which are labeled using an obvious 
notation. The strategy we now follow is to impose that the total scattered field due to 
d{uj) and m(u;) preserves the helicity of the incident field. We first obtain the relationship 
that must hold between d{u) and m{uj) in order for their combined emission to have a well 
defined helicity. Then, we find the conditions that M must meet so that incident fields with 
well defined helicity induce dipoles which produce a scattered field with well defined helicity 
equal to the one of the incident field. 

For the first task, we consider the field emitted by an electric dipole d and a magnetic 
dipole m located at the same point in an infinite homogeneous and isotropic medium with 
electric and magnetic constants (e, /i). We denote by (Ed,Hd) the fields produced by the 
electric dipole d and (E^n, Hm) those produced by the magnetic dipole m. The total fields 
are the sum of the fields radiated by the two dipoles 

E = Ed + E^, H = Hd + H^, (19) 

from which we can obtain the two helicity components: 

G+ = (Ed + E^) + tZ (Hd + H^) , 

(20) 

G_ = (Ed + E^)-zZ(Hd + H^). 
A total field with well defined helicity equal to +1 will have no component of helicity equal 



to -1, thus 

G_ = (Ed + E^) - tZ (Hd + H^) = 0. (21) 

To solve fl2T|) . we use the relations in [ij, chap. 9.3]. According to Jackson, a magnetic 
dipole m produces electric and magnetic fields (E^, Hm) which are related to the electric 
and magnetic fields (Eg, Hg) produced by an auxiliary electric dipole d (note that, for now, 
d and d are not related) in the following way: 



d = -, E, 

c 



-ZHd, Hm = -^^d- 



(22) 



Using ([22]) we turn dZTj) into 



Ed - iZUd = «Ed + ^Hd. 



(23) 



Equation (123|) must be met in all points of space. The solution is hence 

d = -id. 



m 

c 



(24) 



The exact same steps requiring a well defined helicity equal to -1 will result in m/c = d = id. 
We conclude that: 

d = ±1— (25) 

contains the only two cases when an electric and magnetic dipoles at the same point produce 
a field with well defined helicity, respectively equal to ±1. Both types of dipole must be 
present for it. We can now advance to the last part of our program and find the conditions 
on the polarisability tensor M under which the helicity of the incident field is preserved 
in the scattered field due to the induced dipoles in (TT8|) . We proceed by changing our 
representations of the incident fields and the induced dipoles in equation flTSj) in order 
to separate the two helicity components. For the fields, we will use the vector ^ in ([7]). 
For the dipoles, in light of f l25l) . the transformation p± = l/v2 (d ± im/c), separates the 
components that produce fields with well defined helicity. The transformations to obtain Q 
from E and H, and p± from d and m are respectively given by the matrices 



Ti 



1 



V2 



I iZ 
I -iZ 



.To 



1 
V2 



c 



(26) 



We use these matrices to transform equation flTSl) 

d 



To 



m 



TiMT^^Ti 



E 
H 



into 



P^ 
P- 



T2MT{^g. 



(27) 



(28) 



In light ( 128|) . the condition for helicity to be preserved, is that T^MT^"^ must be 3x3 
block diagonal, which then imposes: 

ttdH 



When (129]) is met, we obtain 





"dE = tamH, «mE = 

- - - /i 


Dtain 


P+ 




^dE - i J^^A-R 


P- 




^dE + ^W^^dH 



(29) 



Q. 



(30) 



We conclude that, for scatterers described by their polarisabihty tensors, the relations in 
f l29|) are the necessary and sufficient conditions for helicity preservation, or equivalently, 
duality symmetry. Therefore, and accordingly to our previous arguments, in the context of 
transformation optics and metamaterials equation fl2^ provides a constraint of the electro- 
magnetic response of the meta atoms that does not sacrifice the generality of the achievable 
transformation. Any coordinate transformation can be achieved using only inclusions which 
meet fl29|) . Additionally, the comparison of equations fl30|) and flT^ shows that both e and 
OdE have the same action on both helicity states, while x ^i^nd OdH are the ones responsi- 
ble for helicity dependent transformations. It follows that e can only depend on OdE and 
X can only depend on adH, and, from previous discussions, that adE performs space-only 
transformations and OdH mixes space and time components. 

In this article, we have restricted the polarisabihty tensors of the constituent meta atoms 
in metamaterials used for transformation optics to those that are invariant under duality 
transformations. Duality symmetry, equivalent to helicity preservation, is already a funda- 
mental property of Maxwell's equations in a curved space-time, therefore, the restriction to 
the duality invariant class of meta atoms does not restrict the implementable space-time 
transformations. Additionally, we have shown that the space-only part of the coordinate 
transformation acts equally on both helicity components of the field and depends only on 
the electric-electric term of the polarisabihty tensor, while the part that mixes space and 
time has a helicity dependent effect and depends only on the magnetic-electric term. These 
results can guide the design process between the desired space-time transformation, the cor- 
responding constitutive relations of the transformation medium and its implementation by 
means of metamaterials. Additionally, we have found that for a pair of collocated electric 
(d) and magnetic (m) dipoles to generate a field with well defined helicity equal to ±1 they 
must be related as d = -m or d = — -m, respectively. We have also found the restric- 
tions for helicity preserving (dual symmetric) constitutive relations of inhomogeneous and 
anisotropic media, which include the case of the boundary between two different such media. 
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